Unconventional Multiferroicity in Cupric Oxide 



Pierre Toledano,^ Naemi Leo,^ Dmitry D. Khalyavin,^ Laurent C. 
Chapon,^ Tim Hoffmann,^ Dennis Meier, ^ and Manfred Fiebig^ 

^Laboratory of Physics of Complex Systems, University of Picardie, 33 rue Saint-Leu, 80000 Amiens, France 
^HISKP, University of Bonn, Nussallee 14-^6, 53115 Bonn, Germany 
^ ISIS facility, STFC Rutherford Appleton Laboratory, 
Chilton, Didcot, Oxfordshire, 0X1 1-OQX, United Kingdom 
(Dated: January 25, 2011) 

The magnetic phase transitions reported below 230 K in cupric oxide are analyzed theoretically at 
the macroscopic and microscopic levels. The incommensurate multiferroic and lock-in commensurate 
phases are shown to realize an inverted sequence of symmetry-breaking mechanisms with respect to 
the usual sequence occurring in low temperature multiferroic compounds. The higher temperature 
spin-spiral phase results from coupled order-parameters which decouple at the lock-in transition to 
the commensurate ground state phase. Expressing the order-parameters in function of the magnetic 
spins allows determining the symmetries and magnetic structures of the equilibrium phases and the 
microscopic interactions which give rise to the polarization. 
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In multiferroic compounds such as TbMnOa,^ 
MnW04^ or Ni3V2 08,^ the magnetically induced ferro- 
electric order does not appear directly below the param- 
agnetic phase but across an intermediate antiferromag- 
netic phase via two successive second-order phase transi- 
tions. This typical two-step phase-sequence involves two 
order-parameters associated with the same incommen- 
surate wave- vector. One order-parameter induces the 
paramagnetic-antiferromagnetic transition; it then cou- 
ples to a second order-parameter producing a complex 
spiral-type magnetic order which breaks space inversion 
and gives rise to an electric polarization. Such coupling, 
which expresses the existence of competing magnetic in- 
teractions, tends to lower the transition temperature to 
the multiferroic phase. This phase constitutes in most 
cases the ground state of the system. 

The preceding picture of magnetically driven multifer- 
roics was recently blurred by the discovery of a multi- 
ferroic spin spiral phase in cupric oxide, ^ in which the 
multiferroic and antiferromagnetic phases occurs in an 
inverted order with respect to the standard sequence 
found in low temperature multiferroics. The multiferroic 
phase appears directly below the paramagnetic phase at 
T/V2 = 230K. At Tjsfi = 213K the incommensurate spi- 
ral ordering is replaced by a non-polar commensurate 
spin structure representing the actual ground state of 
the material. The remarkably high transition tempera- 
ture T/V2 has been attributed to large superexchange spin 
interactions, and incidentally related to the property 
of CuO to be a starting material for the synthesis of high- 
Tc superconductors 1^ Here we describe theoretically the 
sequence of transitions in CuO and show that the cor- 
responding transition mechanisms, which both display 
a first-order character, also occur in a reversed order: 
The higher temperature multiferroic phase results from 
the coupling of two order-parameters whereas the low- 
temperature commensurate phase is induced by a single 
order-parameter. Expressing the order-parameter com- 



ponents in terms of spin densities allows to detail the 
magnetic structures of the phases reported from neutron 
diffraction measurements^ri^ and to elucidate the nature 
of the interactions giving rise to the induced polarization. 

The wave- vector kc = (^, 0, — ^) reported for the tran- 
sition to the commensurate phase^ arising at T^i is 
located inside the monoclinic C Brillouin-zone. It is 
invariant by the mirror-plane rriy of the paramagnetic 
space-group C2/cl', corresponding to a two-branch star 

= associated with two bi-dimensional irreducible 
representations (IR's) ri and r2,^^ the matrices of which 
are given in Table |T1 The Landau free-energy express- 
ing the coupling between the corresponding 2-component 
order-parameters {r]i = picosOi,rj2 = pisinOi) and (Ci = 
P2COSO2X2 = p2sin92) reads: 
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Minimization of F truncated at the eighth degree leads 
to eleven possible stable phases for different equilib- 
rium values of pi^p2-,0i and 62- The magnetic sym- 
metries of the phases, numbered from / to X/, are 
listed in Table [III together with the equilibrium values 
of the order-parameters, and the basic translations of 
their fourfold primitive unit-cells. In order to identify 
the order-parameter symmetries associated with the se- 
quence of phases reported in CuO, let us express the 
order-parameter components in function of the magnetic 
spins. Denoting si — the spins associated with the 
eight Cu^+ ions in the magnetic unit-cell shown in Fig. 
[TJa), one can write Si = s^Em + s^bm + 5^c^, where 

am (a,0,c),6rn (0,6, 0), Cm = (a,0, -c), a,6, c be- 
ing the lattice parameters of the paramagnetic C-centred 
unit-cell. Projecting the representation transforming the 
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FIG. 1: (Color online) (a) Different unit cells used 
in the description of CuO. The conventional C-centred 
and primitive paramagnetic unit-cell are represented in 
black and blue solid lines, respectively. The four- 
fold magnetic unit cell is in solid red lines. The 
eight Cu atoms used in our description have the coordi- 
nates 1(1, i, 0), 2(|, |, 0), 3(i, |, i), 4(i, |, -i), 5(-i, |, 0), 
6(-f , 1,0), 7(-l, 1, 1) and 8(-l, 1, -1). (b) Calculated ori- 
entation of the spins in the commensurate magnetic phase of 
CuO (T<213K) along b (red arrows) and parallel to the (a, c) 
plane (green arrows). 



s^{jii = a, 6, c) on ri and r2 one finds three copies of each 
order-parameter (7^1,772) and {C1X2) in function of the 
spin components, which are: 
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Taking into account the equilibrium conditions fulfilled 
by the order-parameter components, indicated in column 
(d) of Table [III the magnetic structures associated with 
the eleven phases listed in Table |II] can be worked out. 




FIG. 2: (Color online) Calculated spin distribution in the 
commensurate approximant of the incommensurate structure 
of CuO (213 K<T<230 K). The sites related by the screw axes 
2i(red and green) of the paramagnetic C2/cl' space group are 
connected by solid lines. 



Using the neutron diffraction data on CuO^^^ then al- 
lows an unambiguous identification of the phases dis- 
closed in this compound. Thus, the commensurate phase 
found below T^i corresponds to phase IV induced by T2 
for = —(^2- It has the magnetic symmetry Pa2i/c 
with the unit-cell basic vectors (2a, 0,0), (0,6,0) and 
(a,0,c) (Fig. m^a)) and the anti-translations Its 
magnetic structure is determined by the equilibrium con- 
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lationships imposed by the ant i- translations, which are 
si -S6, 52 -S5, 53 -54, sj -4. It yields: 
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Fig. mb) shows the corresponding orientations of the 
spins along the h axis and parallel to the (a, c) plane. 
The antiferromagnetic configuration with spins along h 
has been proposed experimentally^^ ^ based on neutron 
diffraction analysis. The model reffects a strong antifer- 
romagnetic super-exchange coupling J along the [101] di- 
rection. In addition, our symmetry-based approach pre- 
dicts that non-zero spin components are allowed in the 
(a, c) plane. The coupling between the in-plane compo- 
nents is, however, ferromagnetic along [101], pointing to 
their probable relativistic origin and therefore too small 
value to be detected in usual diffraction experiments. 

The polar phase VIII of magnetic symmetry Pa2i, 
induced by n + r2, can be taken as a commensurate 
approximant of the incommensurate multiferroic phase 
observed below Tn2- The corresponding equilibrium re- 
lationships, r]f^^ = 7^2^^ and Ci^^ = 



give: 
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TABLE I: Irreducible representations n and T2 of the paramagnetic space group C2/cV associated with the kc — ±(^,0, — ^) 
wave vector star. T is the time reversal operator. 
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Figs. 2 shows the corresponding orientation of the spins 
compatible with the neutron diffraction analysi^ii^ of 
a spiral structure having its envelope parallel to the 
(6*; 0.5a* + 1.5c*) plane in reciprocal space. The an- 
tiparallel orientation of the spins along the Cu-0 chains 
run along the [101] direction, with the strongest antifer- 
romagnetic super-exchange coupling J. Note that using 
phase X of symmetry Pa2 (Table H]) as an approximant 
of the spin-spiral phase of CuO, would lead to the same 
incommensurate magnetic structure which displays the 
point- group symmetry 2yl' . 

Thus, on decreasing temperature below Tn2 the in- 
commensurate phase of CuO results from the coupling 
of two order-parameters which are decoupled at Tati, 
only one of them being associated with the symmetry- 
breaking mechanism giving rise to the low-temperature 



TABLE IL Magnetic space groups (Column (c)) deduced from 
the minimization of the free-energy (Eq. ([I])); Column (a): 
irreducible or reducible representations; Column(b): Stable 
state number; Column (d) Equilibrium values of the order- 
parameters; Column (e) Basic vectors of the fourfold magnetic 
unit-cell. 
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commensurate phase. Figure [3] shows the thermody- 
namic path corresponding to the preceding sequence of 
phases in the theoretical phase diagram of CuO. Note 
that the paramagnetic to incommensurate phase transi- 
tion has necessarily a first- order character since all the 
phases induced by ri + r2 in Table |II] cannot be reached 
directly from the paramagnetic phase across a second- 
order phase transition^i^ 

The polar symmetry of the incommensurate phase al- 
lows an electric polarization component Py, the form 
of which can be deduced from the dielectric free-energy 

2~s ^{ViCi + ^2C2), where 5 is a coupling constant and 

^yy 

the paramagnetic susceptibility. It gives: 
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Eq. (0) is consistent with the linear increase of Py{T) 
on coohng observed below Tim2-'^ It reflects a typical im- 
proper ferroelectric behaviour, which is confirmed by the 



ViPa2/c) 



iC3 



VIIKPJ,) 



Paramagnetic 




C2/cl' 
Tl 



^^^^ (X2 

\lV{PJj/c) 



FIG. 3: (Color online) Theoretical phase diagram associ- 
ated with the free-energy involving the commensurate 
phase and approximant of the incommensurate phase of CuO. 
Solid and hatched curves are first and second-order transition 
curves. Tci-Tc4 are tricritical points. Ti and T2 are triple 
points. The arrow shows the thermodynamic path followed 
on cooling in CuO. 
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upward finite discontinuity of the dielectric susceptibil- 
ity. By contrast the absence of noticeable discontinuity at 
Tn2 for the polarization^ does not allow confirming the 
first-order nature of the transition imposed by symme- 
try and topological requirements.^^ In this respect, addi- 
tional precise experimental studies are necessary to ver- 
ify this important point. Taking into account the order- 
parameter copies Py reads: 



25^ 



and Akz = 



Py = 
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The Sab I Sba and Scb / Sbc coupling invariants contain exclu- 
sively antisymmetric invariants of the s'^s^ — type, 
which represent the Dzyaloshinskii-Moriya (DM) "cy- 
cloidal" interactions coupling the s\ components to the 
or spin-components. The 5ac/Sca terms in Eq. (|6j) ex- 
press symmetric interactions of the 5^5^ + type cor- 
responding to anisotropic exchange between spins, which 
is usullay much weaker than DM terms 

At the phenomenological level the incommensura- 
bility in cupric oxide is related to the existence 
of anti-symmetric (Lifshitz) invariants of the form 

- m'-§t and Ci^ - C2f^(« = x,z) permit- 
ted by the order-parameter symmetries. The free- 
energy of the incommensurate phase reads: $ = 
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of $ yields the incommensurate wave-vector kinc = 



kc + (A/c^,0, A/c^) with Ak^ .~ ^.^_4^^^^ 

Akx{x ^ z), corresponding in CuO to kinc = 
(0.506,0,-0.483). Close below Tn2 the Lifshitz invari- 
ants favour the onset of the incommensurate phase. On 
coohng, the order-parameter magnitude increases and 
fourth-degree terms in Fi influence the stability of the 
system. The first-order lock-in transition to a commen- 
surate phase at Tni coincides with the onset of the 
anisotropic 7^ and Ui invariants in Eq. ([T]) and with a 
cancelling of the Lifshitz invariants. 

In summary, our theoretical analysis of the unusual se- 
quence of phases observed in cupric oxide shows that the 
corresponding transition mechanisms follows in many re- 
spects an unconventional scheme as compared to the situ- 
ation found in low- temperature multiferroic compounds. 
The phases arise in an inverted sequence diCioss first- order 
transitions, the higher-temperature multiferroic phase, 
induced by the coupling of two antiferromagnetic order- 
parameters having a lower magnetic symmetry than the 
low-temperature phase associated with a single order- 
parameter. The order-parameters have been expressed 
in function of the spin variables permitting a theoreti- 
cal determination of the magnetic structures and of the 
magnetic interactions contributing to the electric polar- 
ization. In this respect the polarization appears as a 
purely induced effect essentially due to DM interactions 
between spins in different directions in space. This con- 
clusion, which derives from the improper coupling relat- 
ing Py to the order-parameter components (Eq. ([6])), 
differs from the interpretation, based on DFT calcula- 
tions, that spin canting and polarization mutually stabi- 
lize each other, or from the weak frustration modelii 
proposed for explaining the high transition temperature 
to the multiferroic phase. 
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